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Abstract 

Cosmic string loops are defined by a pair of periodic functions a and b, which 
trace out unit-length closed curves in three-dimensional space. We consider a 
particular class of loops, for which a lies along a line and b lies in the plane 
orthogonal to that line. For this class of cosmic string loops one may give 
a simple analytic expression for the power 7 radiated in gravitational waves. 
We evaluate 7 exactly in closed form for several special cases: (1) b a circle 
traversed M times; (2) b a regular polygon with N sides and interior vertex 
angle tt — 2ttM/N; (3) b an isosceles triangle with semi-angle 0. We prove 
that case (1) with M = 1 is the absolute minimum of 7 within our special 
class of loops, and identify all the stationary points of 7 in this class. 

PACS number(s): 98.80.Cq, 04.30.Db, 11.27.+d 
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I. INTRODUCTION 



Cosmic strings are topological defects that may have formed at phase transitions as the 
universe expanded and cooled . Loops of cosmic string oscillate and emit gravitational 
radiation. This process is of particular importance because most of the observational limits 
on cosmic string networks in the early universe are obtained by considering the effects of 
this gravitational radiation ( and references therein). 

The power emitted in gravitational radiation by a cosmic string loop depends upon 
its shape and velocity. In the center-of-mass frame, a cosmic string loop is specified by 
the position x(t, a) of the string as a function of two variables: time t and a space-like 
parameter a that runs from to L. (The total energy of the loop is fiL where fi is the mass 
per-unit-length of the string). When the gravitational back-reaction is neglected, the string 
loop satisfies equations of motion whose most general solution in the center-of-mass frame 
is 

x(*,<7) = i[a(* + <7)+b(t-<r)]. (1.1) 

Here a(w) = &(u + L) and b(v) = h(v + L) are a pair of periodic functions, satisfying the 
"gauge condition" |a'(w)| = |b(t>)| = 1, where ' denotes differentiation w.r.t. the function's 
argument. The average power radiated by an oscillating string loop is given by 

p = 7 Gfi 2 c, (1.2) 

where G is Newton's constant and c is the speed of light. Quite suprisingly, the dimensionless 
quantity 7 depends only upon the shape of the cosmic string loop determined by the functions 
a and b; it is independent of the length L of the loop From here on, we therefore set 

L = l. 

In a recent paper, we presented a new formula for 7. The formula is an exact analytic 
closed form for any piecewise-linear cosmic string loop 0. Although numerical values of 7 
are easily obtained by that method, the large number of terms makes it difficult to write 
out the analytic formula explicitly for most loops. The conventions and notation used in the 
present work are adopted from this earlier work. In particular, for the rest of this paper we 
use units with c = 1, and use i, j and k to denote a right-handed triad of orthogonal unit 
vectors in flat M 3 . 

In this short paper we evaluate 7 for a particular class of cosmic string loops. For loops 
in this class, the a-loop traces out the following path. When u vanishes, a(tt) also vanishes. 
As u increases, a(-u) moves smoothly up the z-axis with unit velocity k, until u — 1/2. From 
u = l/2 until u = 1, a(w) moves back down to the origin with unit velocity — k. For loops 
in the class that we consider, the b-loop lies entirely in the x — y plane, i.e. in the plane 
perpendicular to the path of the a-loop. See Figure ^ for an example. Note that the b-loops 
in this class are not required to be piecewise-linear. 

The position (|1.1|) of the cosmic string loop in IB? may be easily visualized. Consider a 
string loop at time t = 0. The x and y components of x(cr, 0) are given by (half) the x and 
y components of b(— a). The z component of x(cr, 0) is given by (half) the z component of 
a(cr). Thus, when viewed from large positive z, the x — y projection of the string loop looks 
like the b-loop (scaled by a factor of 1/2). Of course, the string loop does not lie in a plane. 
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As a increases from to 1/2, the z component of x(<r, 0) increases smoothly from to 1/4. 
As a continues to increase from 1/2 to 1, the z component of x(cr, 0) decreases smoothly 
from 1/4 back to 0. As the value of t changes, the z component of a point on the string loop 
oscillates between and 1/4 while the x — y projection of the string loop remains unchanged. 
From this, it can be seen immediately that string loops in this class will self-intersect if and 
only if the b-loop intersects itself. In this paper, self-intersecting string loops are assumed 
not to intercommute, i.e. the strings are assumed to "pass through" one another without 
interaction. 

A short outline of the paper is as follows. In section 2 of the paper, we derive a simple 
analytic formula for 7 for loops in the particular class described above. This formula ex- 
presses 7 as a convolution of the b function with itself. In section 3 this formula is exploited 
to determine 7 for a several particular loops. In the first case, the b-loop is a circle traversed 
M times in the x — y plane. In the second case, the b-loop traces out a regular, N sided 
polygon with interior vertex angle tt — 2itM/N in the x — y plane. In the third case, the 
b-loop is an isosceles triangle in the x — y plane, with semi-angle 9. In section 4 we find 
all of the extrema of 7 for loops in the limited class described above. We prove that the 
minimum value of 7 is attained only by the string loop of the first case with M — 1, for 
which 

r2ir 1 — COS X 

7 = 16 / dx w 39.002454. (1.3) 

Jo x 

This is followed by a short conclusion, in which we also discuss the rate at which the values 
of 7 for the piecewise-linear polygon case with M = 1 approach the limiting value ( |1.3| ) of 
7 for the perfectly circular case. 



II. FORMULA FOR THE 2-PLANE CASE 

In this section we find a simple general formula for 7, valid for our particular class of 
cosmic string loops. The string loops in this class are defined by a- and b-loops with the 
following properties. The a-loop is composed of two straight segments joined at "kinks" 
where the tangent vector &'(u) is discontinuous. This loop is taken to lie along the z-axis 
with one kink positioned at the origin. (The parameter u is chosen to be zero at this kink.) 
The other kink (at u = 1/2) is positioned above the first kink and has coordinates (0,0,1/2). 
The b-loop is constrained to lie in the plane perpendicular to the a-loop (i.e., the x — y 
plane). The formula for 7 found in this section will be valid for all cosmic string loops in 
this limited class. 

The starting point for the calculation in this section is the following equation, which is 
derived in section 3 of reference 0, 

7 = 4/ du dv du / dv^{u,v,u,v)D e{At)5{{At) 2 - |Ax| 2 ) . (2.1) 
Jo Jo Jo J-00 1 J 

In this equation, 5 is the Dirac delta function and e(x) = 26(x) — l where 9{x) is the unit step 
function. The u and v integrations are over the entire world-tube swept out in space-time by 
the string loop as it oscillates. The u and v integrations are over a region on the world-tube 
swept out by the string loop during a single oscillation. The functions At = (u+v — u — v)/2 
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and Ax = (a(u) + b(t>) — a(tt) — b(v))/2 describe the temporal and spatial separation of 
the two points on the string world-tube with coordinates (u, v) and (u, v) respectively. The 
function if) is defined by 

i/;(u,v,u,v) = -[(a'(u)-a'(u) - l)(b'(v)-b'{v) - 1) + 

( a!(u)-b'(v) - l)(b»-a'(u) - 1) - (a!(u)-b'(v) - l)(b'(v)-a'(u) - 1)], (2.2) 

and D is the linear differential operator 

D = U(u, v, u, v)d u + V(u, v, u, v)d v — U (u, v, u, v)du — V(u, v, u, v)dv- (2.3) 

The functions U, V, U and V are determined by requiring that D satisfy the four equations 

DAt(u,v,u,v) = 1, DAx(m, v, u, v) = 0, (2.4) 

which may be written as the matrix equation 

/ 1 1 1 1 

K a!{u) b'(v) a'(u) b»'(t5) 



V 

u 



(2.5) 



The generic solution to this set of equations is given in equation (3.12) of 0. 

We restrict our attention to the limited class of loops described above for the remainder 
of the paper. With this restriction, ( |2.1| ) will lead directly to a simple formula for 7. The 
function a(u) for this limited class of string loops may be written as 



a(it) 



uk for u G [0, i), 

(1 - u)k for u e [~, 1), 



(2.6) 



where k is the unit vector along the z-axis. From this equation and the fact that the b-loop 
lies entirely in the plane perpendicular to k, it follows immediately that ip takes the form 



ip(u, v, u, v) 



for u e [0, |), u E [0, |) or u G [§, 1), u G [|, 1) 

|(1 - b'(v)-b'(v)) for u G [0, |), tZ G [|, 1) or it G [|, 1), t2 G [0, |). 



(2.7) 



In addition, for u G [0, |) and u G [|, 1), equation (|2.5| ) for the coefficients in the differential 
operator D becomes 



/l 











V 






u 











(2.8) 



with the obvious solution U — U — 1, V — V — 0, i.e, 



4 



D = d u - du. 



(2.9) 



It is clear that the same operator also satisfies equation (|2.4j ) for u G [|, 1) and £t G [0, |). 
(It is interesting that this solution works whether or not the determinant of the matrix in 
Q2.5|) vanishes.) Substituting ( [2.7|) and (|2.9j ) for ip and Z) into Q2.1|) , we find that 
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cfc / dv (l-b'{v)-b'{v)) x 



/ du + du du \{d u -du) e{At)5{{Atf - |Ax| 



(2.10) 



For the specific form of a(w), the second integral in curly braces immediately reduces to the 
first under the change of variables u — 1 — u', u — 1 — v! . In addition, we may perform 
the same change of variables in the second term of the derivative operator. Doing so leads 
directly to 
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where 



Q{v,v) 



dv / dv {l-b'(v)-b'{v))Q(v,v), 



e(At)8((At) 2 - |Ax| 



du / du <9„ 



(2.11) 



(2.12) 



The next step will be to carry out the integrations over u and it in Q{y,v). 

The partial derivative in ( p.!2|) allows the integral over u to be performed immediately. 
In doing so, the function in square brackets in ( p,12j ) must be evaluated at both the upper 
and lower limits of the u integration, yielding 



Q{v,v) 



du 



iAt)5((Atf 



|Ax| 



2 

u=0 



T(v,v) - B(v,v), 



(2.13) 



where T is the integral of the quantity in square brackets evaluated at the top limit, and B 
is the integral of the quantity in square brackets evaluated at the bottom limit. Consider 
first the top limit, T(v,v). Because u G [|, 1) we have 



1 _ 
At = - — u + v — v, 

Ax = (u - -)k + b(v ) - b(v). 



Introducing u' — u — |, the integral over u in T{v , v) can be written as 



1 r 



T(v,v) = 

Turning to the lower limit, B(v,v), we have 



du' e(v-v- u')8 ( ^ [-2u'(v -v) + (v- vf - (b(v) - b(d)) 2 



At = —u + v — v , 
Ax = (u - l)k + b(v) - b(v). 



(2.14) 
(2.15) 



(2.16) 



(2.17) 
(2.18) 
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We now use the freedom to rewrite (|2.11|) as an integral with v running from -oo to oo 
and v running from to 1. (This may readily be seen as follows. Rewrite the integral 
over v as a sum over n from — oo to oo of an integral from to 1 of the integrand with v 
replaced by v + n. Then use the explicit form of the integrand and the periodicity of b to 
combine the n with the v to transform the v integral into an integral from — oo to oo. J7|) 
If we now make the simultaneous changes of variable v — v ', v — v' + 1 and u = 1 — u' we 
recover exactly fl2.11| ) (with v and v replaced everywhere by v' and v') except that the sign 
of the argument of the e-function has changed. However, since the e-function is odd and we 
are dealing with the bottom limit in ( |2.13| ), we recover exactly ( 2.1tj| ) for —B(v,v). Hence, 
Q(v,v) = 2T(v,v). 



Using (|2.13|) , (12.16 ), and T = —B, Q(v, v) may be rewritten in the form 



(v - vf - (b(v) - b(v)Y , 
Q(v,v) = -. du' e{v - v - u')8 [u'-£ — ~r~ ~z — - 1 . (2.19) 



v\ Jo \ 2{y — v) 

This integral may be carried out immediately, yielding 

QM = _* Uv-WHm-m)') x 

\v — v\ \ 2(v — v ) J 

J (v - vf - (b(v) - b(D)) 2 \ (\ [v-vf- (b(v) - b(v)) 2 \ 
\ 2{v-v) j \2 2{v-v) )' 1 ' ' 

This answer may be considerably simplified. 

First, because the numerator of the argument of the e-function in Q(v,v) is manifestly 
positive, the function is just equal to e(v — v). This term simply has the effect of removing 
the absolute value signs from the denominator of the prefactor. Next, consider the two 
^-functions appearing in ( |2.20|) : 



(1) \v — v'\ is the arc-length along the b-loop from b(v) to b(v) (traversed in a particular 
sense possibly more than once) which must be greater than the straight line distance \h(v) — 
b(v)\ between the two points. Thus again the numerator is a positive function and the first 
^-function is just equal to 6{v — v). 

(2) The argument of the second ^-function can be written as 

(l-(v-v))(v-v) + (b(v)-b(v)) 2 ^ 
2(v — v) 

We may neglect the denominator which is always positive when the first ^-function is non- 
vanishing. In addition, we need only consider the behavior of the numerator for v — v > 0. 
First we note that for (v — v) G (0, 1) it is manifestly positive. On the other hand for 
(v — v) > 1, we note that (v — v) is the arc-length along the b-loop from b(v) to b(v) 
traversed in a particular sense (possibly more than once) while (v — v) — 1 is the arc-length 
along the b-loop from b(v) to b(v) traversed in the opposite sense (possibly more than 
once). Thus (v — v) > \b(v) — b(v)\ and (v — v) — 1 > \b(v) — b(v)\ and so 

(1- (v -v))(v -v) + (b(v) -b(v)) 2 <0 for (v - v) > 1. (2.22) 

Combining these results we see that ( [2.20|) simplifies to yield 
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Q(v,v) = -^—6{v- v)B{l-{v -«)). (2.23) 
[V — v) 



Replacing Q(v,v) in ( |2.11| ) with ( j2.23j ), we obtain our principal result, 



7 = 16 / dv / dv 

JO Jv-l 



l-b'(v)-b'(v) 



[V — V) 



= 16 fdv f 1 fc l-b»-b>-,) 

Jo Jo X 

In the next section we shall evaluate this remarkably simple expression for a number of 
particular loops. 



III. RESULTS 
A. Circular b-loop 

Our first application of formula (|2.24j) for 7 is the case where the b-loop winds M times 
around a perfect circle (see Figure [I]). Here M is any positive integer. In this case, the 
function b(t>) may be written as 

b(v) = — -r[cos(27rMw)i + wa.(2irMv)j\. (3.1) 
2ttM 



The dot product appearing in ( 2.24 ) gives 



b'(v) ■ b'(v -x) = [cos(27rMw) cos(2vrM(w - x)) + sm(2irMv) sin(2vrM(t; - x))} 

= cos{2nMx). (3.2) 

Thus, 7 is given by 

f 1 , f 1 , 1-cos(2ttMx) _ f 1 , 1-cos(2ttMx) 
7 = 16 / au / ax = 16 / ax 



Jo x jo x 

p'2tvM 1 COS T* 

= 16 / cix = 16 [C + ln(27rM) - Ci(27rM)] (3.3) 

Jo X 

where C = 0.577216 ... is Euler's constant and Ci(x) is the cosine integral function defined by 
equation (5.2.2) of reference ||. The numerical value of 7 for this string loop is 7 ~ 39.002454 



when M = 1. In section [TV] we show that this is the minimum value of 7 for any loop in the 



limited class of loops for which the formula ( p. 24 ) is valid. 



B. Polygon b-loop 



As the next application of formula ( |2.24| ), we consider the case where the b-loop takes 
the shape of a regular, iV sided polygon (see Figure ||). Because the length of a given side 
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of the polygon is l/N, it is convenient to break the integrations in ( [2.24 ) into iV equally 
spaced intervals, 

7 = 16 E * IT'" ^- b ' {v) : b ' (V - X) . (3.4) 

i,j=0 



i/JV J j/N x 



For each i and j, the t> and x integrations are over specific segments on the polygon. The 
integrations may be simplified by introducing a new pair of coordinates, v and x , which 
will each be in the range [0, l/N] : 

u = «--!, x = x-jj. (3.5) 

With this change of variables, 7 becomes 

7 = 16 V / V " A <a l-^ + ")-b'(^ + «-i/iV-x) - (3 

,4^nio io i/N + x v ; 



ij=0 



The vector b'(z/iV + 7;) is a constant vector (tangent to the zth segment on the polygon) for 
each value of i. When x < v, the vector h'(i/N + v — j/N — x) will be the constant vector 
tangent to the (i — j)th segment on the polygon. When x > v, it will be tangent to the 
(i — j — l)th segment. Thus, the x integration may be broken into two ranges where the dot 
product h'(i/N + v) ■ h'(i/N + v — j/N — x) is constant in each range. 

V* f r >- 1 - b'(z'/iV + • b'(z/iV + v - j/N - x) 



7 = 16 / \ I dx 



■i/jy ^ 1 _ b '(ViV + g) ■ b'(z/iV +^- j/jV - x) . 
' j/N + x 1 1 K,! 



i,j=0 



For each value of i and j, the constant dot products in (|3.7|) may be evaluated. Because 
h'(i/N + v) is tangent to the zth segment on the polygon and h'(i/N + v — j/N — x) is 
tangent to either the (i — j)th segment (for v > x) or the (i — j — l)th segment (for x > tJ), 
it follows that their dot product is given by 

{cos(^j) for v > x, 
(3.8) 
cos (^r(j - !)) lorv<x, 

where 2ir/N is the angle through which one must rotate a segment of the polygon to bring 
it parallel to the next segment. The r.h.s. of (|3.8j ) is independent of i. Thus, the sum over 
i simply gives a factor of N. The integrations over x may be done giving 

7 = 16iV£ / dv\(l-cos(—j))(\n(j/N + v)-\n(j/N) / 

+ (l - cos(^(j + 1))) ( ln(j/iV + l/N) - \n(j/N + «))}. (3.9) 
The v integration may now be carried out. After combining terms, this yields 
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TV— 1 r O O I 1 

7 = 16 E (l + 3 cos(-(j + 1)) - (J + 1) C ^(-J)) In (^-) 



.27T /27T 

+ cos(— j) - cos(— (j + 1)) 



+ 16(l-cos(^)). 



(3.10) 



Finally, we note that all the terms not multiplied by In (^y=) in ( |3.10| ) cancel exactly in the 
sum over j. Thus, we arrive at the final form 



N-i 2ir 
7 = 16 £ (l + 3 cos( — (j + 1)) - (j + 1) cos(— j)) In ( 



(3.11) 



or equivalently 



- = 32 ( 1 - cos(-)] (-Win JV + E J «»(— / 



(3.12) 



For the first few values of iV, equations ( |3.11| ) and ( |3.12| ) give 

' 64 In 2 for N = 2, 

7=S 721n3-481n2 for N = 3, 

641n2 for iV = 4. 



(3.13) 



Note that the N = 2 case is identical to the Garfinkle and Vachaspati case (equation (3.9) 
of reference 0) with a = tt/2, and has the same value of 7 as they obtained. Numerical 
values for 7 are shown as a function of N in Figure [| 

It is interesting to note that ( |3.11| ) may be trivially modified to obtain a formula for 7 for 
a set of (self-intersecting) loops related to the polygon case. Consider again a b-loop which 
is composed of N straight, equal length segments. However, in this case, instead of placing 
each segment at an angle of 2ir/N relative to the previous segment (and thereby getting an 
N sided polygon), each successive segment is placed at an angle of 2Mtt/N relative to the 
previous segment (where M is a positive integer not equal to N). This causes the b-loop to 
wind M times around the origin. The only change required in ( |3.11| ) is to replace the factor 
2tt/N appearing in the cosine terms by 2M7r/N, 



n-i 



^ 1 ,2nM , XN . , ,2-kM .x, /? + l 

7 = 16 E (1 + j cos(— (j + 1)) - (j + 1) cos(— j)) In (^— 



(3.14) 



or equivalently 



7 = 32 1 - cos( 



,2vrM 



\/l N ^ 
)(-N\nN+ cos 

' \* j=2 



,2-kM 
- N " 



(3.15) 



If M = 1, we have the polygon case considered above. If for instance, M = 2 and N — 5, 
the b-loop would have the shape of a pentagram (see Figure [|). 
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C. Isosceles b-loop 



Our final application of ( J2.24|) is the case where the b-loop takes the shape of an isosceles 
triangle. This is a one-parameter family of loops where the parameter is the half-angle 9 
between the two equal length sides of the triangle (see Figure |5]). The length of the two 
equal sides of the triangle is denoted by l\. The length of the base of the triangle is denoted 
l 2 - Using the constraint 2li + l 2 = 1, one can easily write l\ and l 2 in terms of the semi-angle 
9 of the triangle, 



1 



2(sin0 + 1)' 



l 2 = 21 1 sin 9. 



(3.16) 



The parameter along the b-loop is taken to run from to l\ along the first side of the 
triangle, from li to 2l\ along the second side, and from 2l\ to 1 along the third. If we label 
the unit tangent vectors to the first, second and third sides of the triangle by b' 1; b' 2 and b 3 
respectively, the dot products between the various sides are given in terms of 9 by 



b[b' 2 = -cos(20), 
b;.b'3 = b' 2 -b^ = - 



sin 9. 



(3.17) 



The dot product of any unit tangent vector with itself is simply 1. 

We are now in a position to use equation (|2.24j) to find 7 for this family of string loops. 
The integrations over v and v may be broken into ranges corresponding to the three sides 
of the triangle, yielding 



7 = 16< 



dv 



h-h I - 
dv — 

v-l V — V 



+ 



2/1 



dv 



dv 



2/1 



bi • K 



dv — 



b[ ■ b' 2 



+ 



h-h 

k ■ b: 3 



v — V 



dv- 

-h v -v 



+ 



1 — h' • b' 

dv 1 1 



v —v 



l * , l-b'-b' 
dv 

v-l V — V 



+ 



1 — h' • b' 

dv 2 3 
-h v — v 



fh 1 - b' 2 • b; p> 
+ / dv 2 - r - J - + / dv 

Jo v — V Jh 



1 - b 2 • b 2 

v — V 



dv — 
v-l v — v 



b' 3 • b 3 



■'i , l-b'-b' 
+ / dv z - l - 

v — V 



,„ 1 - b' 3 • b' 2 

+ / dv 

Jh v — v 



+ I dv- 

2h 



-b 3 b 3 

v — V 



(3.18) 



The numerator of each integral in ( |3.18| ) has the form 1 — • b^- . Each of these numerators 
is now constant. If i = j, the numerator (and hence the integral) vanishes. When % 7^ j, 
equation fl3.17|) may be used to write the dot products in terms of 9. The integrals in fl3.18| ) 
are now easy to evaluate and yield 



7 = -16{ ( 1 + cos(20))(2(Zi + k) Hh + h) ~ h Hh) 
+ (l + sin0)(2/ 1 ln(2/ 1 ) + / 2 \n(l 2 ) + l 2 ln(/ 2 ) - 



f 2/iln(/i) -2/i ln(2/i)) 
2/ 1 ln(2/ 1 ))}. (3.19) 
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Using ( |3.16| ) to write l± and li in terms of 9 and simplifying gives the final form 
7 = 32 [ 2 In 2 cos 2 9 - sin 2 (#) m(sin 0) + (2 - sin 5) (1 + sin 9) ln(l + sin 9) 



1 - sin 0) (1 + 2 sin 9) ln(l + 2 sin I 



(3.20) 



This formula for 7 is plotted in Figure |6] for a range of angles 9. It should be noted that 
( |3.2U| ) correctly reduces to the M — 1, N = 2 result ( ^.14[ ) of the previous section when 
= 0° or 90°, and to the M = 1, iV = 3 result when = 30°. 



IV. MINIMUM 7 LOOP AND STATIONARY POINTS OF 7 



In this section, we obtain the minimum value and stationary points of 7 for all loops in 
our limited class. Since b(t>) is a periodic function on the interval v G [0, 1), its derivative 
w.r.t. v is also a periodic function, and may be expressed as the Fourier series 

b»= £ e 2mnu (a n i + b n i). (4.1) 

n=— oo ^ ' 

Since the function h'(v ) is real, one has a ra = a^ n and b n = b*_ n with * denoting complex 
conjugation. The periodicity of b also implies that a® and bo vanish, since 

1 b'(v)dv = b(l) - b(0) = a l + b j = 0. (4.2) 

For the remainder of this section, we use the notation J2' n to indicate the sum from minus 
infinity to infinity with the n = term excluded. 

The Fourier coefficients a n and b n are not arbitrary; they satisfy an infinite number of 
constraints which follow from the gauge condition that the length of the vector |b'(u)| = 1. 
One of these constraints may be expressed in a useful "integrated" form: 



\b'(v)\ 2 dv 



EE e 

n m 



2iri(n+m)v 



b n b m )dv — E ( a n«n + b n b* r 



(4.3) 



where in the final line we have used the orthogonality of the exponential functions on the 
unit interval and the reality conditions a n = a*_ n and b n = b*_ n . One may define a non- 
negative real quantity c 2 = 2(a„a* + b n b* n ). Note that for any cosmic string loop one has 



,.2 

-n 



cz = c „ , and that the normalization condition 



E< 

71=1 



is implied by the integrated constraint (|4.3|) . 

Now consider the value of 7 given by equation (|2.24| ), 



7 



16 f 1 dv f 1 — [1 -b'(v) -b'(v-x)] 
Jo Jo X 



16 



1 dx 



x 



dv 



E' E' e 2m « n+m >~ m ^ (a n a m + b n b T , 



y 1 dx 
16 / — 

x 



E'' 



> 2 ™ x (a n a* n + b n b* n ) 



(4.4) 



(4.5) 
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For any cosmic string loop, the function b satisfies the integrated constraint ( f4.3|) . Substi- 
tuting this for the first term in the previous equation ( |4.5|) , one obtains 

„\ /i „2ninx\ co 

7 = 16 Y! dx K - L {a n a* n + b n b* n ) = £ X n c 2 n . (4.6) 

r, JO X n = 1 



This is a quadratic form for 7; its "eigenvalues" are 



A„ = 163? 



1 (1 -e 



2mnx\ 



dx 



2-ku (I _ cosx ) 



16 / dx 

x 



>0 X 

= 16[C + ln(27rn) - Ci(27rn)] 

roc te~ 27Tnt 

= 16[C + ln(27rn) - / -dt } 

Jo 1 + 1 2 

« 1Q[C + \n(2nn) - — ^ + 0{\) }, (4.7) 

where C = 0.577216... is Euler's constant and Ci(x) is the cosine integral function defined 
by equation (5.2.2) of reference ||. From the integral forms it is clear that the eigenvalues 
increase without bound: A n < A n+ i. The first few eigenvalues are approximately 

(Ai, A 2 , A 3 , • • ■) « (39.0025, 49.8297, 56.2636, 60.8473, 64.4086, 67.3208, ...). (4.8) 



From the quadratic form of 7 (4.6) one may easily find the absolute minimum and stationary 
points of 7. 

Since any cosmic string loop in our limited class has J2n°=i c n = 1 an d 7 = S^Li c ^A n , it 
is immediately clear that the minimum value of 7 is attained when cf = 1 and all other cjL^ 
vanish. (The problem is identical to that of minimizing the energy of a quantum mechanical 
system with ortho normal eigenstates \n >, where the general normalized state is denoted 
J2n c n\ n > and the energy of the n'th state is A n .) The loop with c\ = 1 is precisely the 
M = 1 case considered in Section 3A, where the b-loop is a circle traversed once as v 
increases from zero to one. The remaining stationary points of 7 are immediately seen to be 
loops with c\ = 1 for some integer k and all other coefficients vanishing. This corresponds 
to the M = k case considered in Section 3A where the b-loop is a circle traversed k times 
as v increases from zero to one. 



V. CONCLUSION 

In this paper we have found a simple formula for 7, the power radiated in gravitational 
waves, for a limited class of cosmic string loops. We have applied this formula to a variety 
of loop configurations within the limited class and have shown that the minimum value for 
7 in this class is attained only when the b-loop has the shape of a perfect circle (traversed 
once) . 

The results for 7 in the piecewise-linear cases (polygon and isosceles b-loops) have been 
checked against the general analytic formula given in reference , using a symbolic manip- 
ulator ( Mathematica ) to evaluate the formula of reference M. We have also used these 
exact results to test a floating-point implementation (written in C) of the exact formula of 
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reference 0. This floating-point implementation was described in some detail in [7]] and 
is publicly available. The tests have shown that the floating-point implementation typi- 
cally yields results accurate to between 6 and 8 decimal places, the difference being due to 
accumulated round-off error. 

In reference |7j we presented an analytic formula yielding 7 for piecewise-linear loops, 
and argued that any smooth loop could be well-approximated by an iV-segment piecewise- 
linear loop with N sufficiently large. The exact formula in this paper allows us to determine 
the rate of convergence for large N, at least for a limited class of loops. The results for the 
circular b-loop case considered in section |TTT| A (with M — 1) may be used to examine how 
well smooth loops are approximated by piecewise-linear loops. The regular polygon b-loops 
(M = 1) of section are increasingly circular in shape as N — > 00. The difference between 
the 7 value for a cosmic string loop with a b-loop in the shape of an N sided polygon (7^) 
and the result from the case where the b-loop is a perfect circle (700) is shown as a function 
of TV in Figure 0. The piecewise-linear result quickly converges to the smooth result (see 
also Figure |3|). The convergence goes like 1/iV 2 for N > 10. 
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FIGURES 



FIG. 1. A 2-segment a- loop perpendicular to a circular b-loop with M = 1. The cosmic string 
loop defined by these a- and b-loops has 7 ~ 39.002454. This is the minimum value of 7 attained 
by any loop in the limited class which we consider. 

FIG. 2. A 2-segment a-loop perpendicular to a b-loop in the shape of a regular polygon with 
M = 1 and N = 5. The cosmic string loop defined by these a- and b-loops has 7 ~ 43.100206. 

FIG. 3. Values of 7 for cosmic string loops with 2-segment a-loops perpendicular to b-loops in 
the shape of regular N sided polygons with M = 1. These 7 values converge like 1/N 2 to the value 
of 7 for the circular b-loop case. 

FIG. 4. A 2-segment a-loop perpendicular to a b-loop in the shape of a regular polygon with 
M = 2 and N = 5. The cosmic string loop defined by these a- and b-loops has 7 ~ 56.980278. 

FIG. 5. A 2-segment a-loop perpendicular to a b-loop in the shape of an isosceles triangle with 
semi-angle 9. The semi-angle 6 lies in the range 9 G [0, vr/2]. 

FIG. 6. Values of 7 plotted as a function of the semi-angle 6 for the cosmic string loops 
described in Figure 5. The cases 9 = and 9 = 90° yield the same result, 7 = 64 In 2 ~ 44.361419, 
as Garfinkle and Vachaspati equation (3.9) of reference || with a = ir/2. 

FIG. 7. The convergence of A = 77V — 7oo, where 77V is the power radiated by a cosmic string 
loop with a 2-segment a-loop perpendicular to an N sided regular polygon (M = 1) b-loop, and 
7oo is the power radiated when the b-loop is a circle (M = 1). The value of A converges to zero 
like l/N 2 for large N. 
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